Binding of charmonium with two- and three-body nuclei 



O 
O 

(N 

C 



i 



=5 



X 



V.B. Belyaev, N.V. Shevchenko* 
Joint Institute for Nuclear Research, Dubna, 141980, Russia 

A. Fix 

Institut fur Kernphysik, Johannes Gutenberg- Universitdt Mainz, D-55099 Mainz, Germany 

W. Sandhas 

Phystkahsches Institut, Universitdt Bonn, D-53115 Bonn, Germany 
(Dated: February 9, 2008) 

The energies of the rj c d and 7^ c 3 He bound states are calculated on the basis of exact three- and 
four-body AGS equations. For the r/ c N interaction a Yukawa-type potential has been adopted. The 
calculations are done for a certain range of its strength parameter. The results obtained are quite 
different from calculations based on the folding model. 
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I. INTRODUCTION 



The behavior of mesons in a nuclear medium has attracted much attention during the last years [|. We mention 
the existence and properties of deeply bound meson-nuclear states Q as well as the possible restoration of chiral 
£S| ■ symmetry in the nuclear environment |3j , a point closely related to modifications of widths and masses of the mesons 
' propagating through nuclear matter |4j . 

It should be emphasized that all these phenomena arc usually described in terms of effective degrees of freedom 
including mesons, nuclcons, and isobars. The manifestation of more fundamental quark-gluon degrees of freedom at 
low energies is very unlikely because of the strong mismatch between the nuclear energy scale (few tenth of MeV) and 
«) ■ the QCD energy scale (few hundreds or thousands of MeV) . 

However, as pointed out in [5j with regard to charmonium-like mesons, for example T) c , which do not contain u- and 
• e?-quarks, the main contribution to the meson-nucleon interaction is expected to originate not from effective, but from 
fundamental QCD degrees of freedom, namely from the few-gluon exchange. In this way, using the Pomeron exchange 
model, one can deduce a Yukawa-type potential for the r^-nudeon system Following these ideas, the authors of 
were able to explain the pp-spin correlation at the threshold energy of charm production. Such an analysis allows 
one to restrict the uncertainties of the parameters of the charmonium-nucleon potential. Independent theoretical 



indication on the existence of attraction between charmonium and nucleons has been obtained in [7j . In case of an 
infinite-mass limit of the heavy quark, these authors have obtained a binding energy of J/ip in nuclear matter of 
the order of 10 MeV. For the (2S) charmonium state this approach predicts an even stronger binding. Despite the 
above indications, the data of the charmonium- TV interaction are still rather scarce. Bearing this in mind, it appears 
reasonable to look for possible formations of the corresponding mesic nuclei, similar to the 77-nuclear systems which 



were extensively discussed in the literature (see e.g. |8|] and references therein). In particular the study of the reactions 

p + A ->p + p + A' -> T) c + A' (1) 

has been planned within the PANDA project The production of r] c in electromagnetic processes has already 
been started in the CLIO collaboration [l0|. As explained in p|, the precise knowledge of the binding energy of 
meson-nuclear systems is of crucial importance for their experimental observation. 

First attempts to calculate the binding of ?7 c -mcsons with light nuclei have been made in Ref. Q on the basis of 
variational calculations, and in Ref. 11] where the folding model was used. Both calculations involve a Yukawa form 
for the basic rj c N interaction. Our results for the binding energies of rj c d and r\ c 3 He, presented below, are obtained 
within the AGS theory 0,0]. Similar to 0,^3 we use the Yukawa type rj c N potential 

K )c Ar(r) = -a— (2) 
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with a = 0.6 GeV. The parameter a is varied between 0.4 and 0.6. Since the two-body potentials act essentially 
via their s-waves contributions, we take into account only s-wave parts of the corresponding rj c N and NN scattering 
amplitudes entering the kernel of the AGS equations. 



II. FORMALISM 



Let us start by giving a brief description of the 3-body rj c NN formalism. Its main ingredient is a separable 
representation of the driving rj c N and NN potentials 

N 

V a =J2\xi)(xi\, ae{ Vc N,NN}. (3) 

i,3 

To outline the important aspects of this approach, we use for the moment an oversimplified representation of the rj c N 
potential keeping only the first term in the sum J2Jl , while the actual calculations are performed with a more realistic 
representation containing six terms (we will discuss the sensitivity of our results to the number of terms kept in the 
separable expansion). For the s-wave triplet state 3 Si of two nucleons we employ the separable rank-one version of 
the Paris potential of Ref. • 

The AGS equations for the rj c NN system read jg 

X u (z) = 2Z 21 (z)r 2 (z- ^)X 21 (z) (4) 

2 2 

X 2 i(z) = Z2i+Z2i(z)n{z--^)X n {z) + Z2 3 (z)r2{z--^-)X2i(z). 

L\X 2 

The functions r a (z), Z a p, and X a p (a, (3—1,2) can be expressed in terms of the AGS transition operators Up a and 
form factors |x Q ) appearing in J2J via 

Zp a (z) = (1 - 6 a p){xp\Ga(z)\xa) (5) 
Xp a (z) = {xp\G (z)Up a G (z)\ Xa ) (6) 

where Go(z) is the free rj c NN propagator. As mentioned, for the NN subsystem we use a separable representation 
of the Paris potential 

6 £, 

V N N(k,k') = -g(k)g(k') , where g(k) = ^ fc2 ' , (7) 

i—l 2 

with parameters C; and (3i listed for the 3 Si configuration in [T^ . 

For the 7? e iV interaction we first consider the rank-one representation of the r\ c N potential J2J within the Bateman 
method The corresponding general expression reads 

JV 

V^ N] (k, k') = V Q (k, Si ) d"- 1 V Q (sj,k'), (8) 

with dij = Vo(si, Sj). Here Si is an expansion parameter and 

1 f°° 

V (k,k') = — / ]() {kr)V{r)j Q {k'r)r 2 dr (9) 

the s-wave Fourier transform of the potential V(r). Using expression J5J) for the potential V(r) in the integral © 
and taking N — 1 and si = in the expansion (jSJ) one can readily obtain 

If the Bateman method is applied to a negative potential, the next terms in the separable expansion tend to increase 
the attraction. Therefore, the neglect of the higher order terms in the sum (jHJ will result in underestimation of the 
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FIG. 1: Binding energy of r\ c d system, as function of the strength a of the Yukawa potential lp)l. The results obtained with the 
Bateman and Hilbert- Schmidt representation of the r/ c N potential are presented by the dashed and the solid line, respectively. 



absolute value of the rj c d binding energy. An obvious advantage of the Bateman method is that it provides analytical 
expressions for all terms in (jSJ. 

It is worthwhile to note that the potential (|10|) predicts an antibound (virtual) state in the rj c N system lying at 




= -6.1MeV, (11) 

where /i is rj c N reduced mass. The corresponding scattering length 

a VcN = 4^— « 0.97 fin (12) 
— — a 

is a bit larger than the r/N scattering length predicted by the modern analyses (see e.g. compilation in 0])- Fur- 
thermore, there is no absorptive part in the r] c N potential, which in the r)N case effectively increases the repulsive 
properties of the interaction. The most important difference between rjN and r] c N is the larger mass of the r\ c meson, 
which should increase the role of attraction in the rj c few-nucleon dynamics. According to this, one can expect that 
the rj c meson should be more tightly bound to the nucleons than the r\ meson. The energy of the rj c d bound state, 
calculated with the potential JTU)> as a function of the strength parameter a, is presented in Fig. 1 by the dashed line. 

As a next step we employ the Hilbert-Schmidt method for the separable representation of the potential (J2J. The 
corresponding expansion reads 

iv 



Vr, eN (k, k') = Kv n {k)v n {k') , (13) 

n=l 

where the form factors v n (k) and strength parameters A„ obey the homogeneous integral equation 

v n (k) = i- Jvr,Mk,k') ^% k ' 2 2 ^ • (14) 



B = was chosen in the actual calculation. 



III. RESULTS 



The r) c d binding energy calculated for the potential (|13fl with N = 6 terms is presented in Fig. 1. As expected, the 
last model predicts stronger binding than the leading term in the Bateman expansion (jHJ. Furthermore as one can 
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see, in contrast to the results of i and [12, 

we find only one bound state lying in the region —2 MeV < Ef, < —5 
MeV, when the parameter a is varied between 0.5 and 0.6 GeV. 



TABLE I: Dependence of the r/ c d binding energy Eb on the number of separable terms N in the expansion 1131 . Parameter a 
is equal to 0.55. 



N 


1 


2 


3 


6 


E b [MeV] 


-3.429 


-3.251 


-3.217 


-3.196 



The results in Table I show the accuracy of the Hilbert-Schmidt method for the rj c d binding energy. Evidently, 
already the first two terms in the sum 113(1 provide quite a satisfactory approximation. Going up to N — 6 results in 
only 1.7% increase of \Eb\. 

Now we turn to the four-body r\ c 3 He system. For the calculation we use the formalism developed in Ref. 01 
and Ref. |18| for the interaction of isoscalar-pseudoscalar mesons with three-body nuclei. Since the formal part of 
the problem is quite involved we do not present it here and refer the reader to these references. It should only be 
recalled that the calculations are based on separable representations of the two- and three-body kernels of the basic 
AGS equations so that, as a consequence, the resulting integral equations have the same structure as the three-body 
equations Q). 





TABLE II: 


The binding energy of r\ c 3 He system (in MeV). 




a 


4-Body calculations 


Variational calculation[5j 


Folding model ^Tj 


0.4 


-1.3 


-3.0 




0.6 


-14.5 


-19.0 


-0.8 



Let us compare our results for the binding energy of rj c 3 He system with variational and folding model calculations 
|ll|. One should emphasize that the variational method used in [j| is not an ab initio many body treatment. Rather 
it is a model of meson-nuclear interaction, postulating a Yukawa form for the meson-nuclear potential with a range 
depending on the number of nucleons in the nucleus. This model, as well as the folding model, takes into account the 
important effect of smearing the meson-nucleon interaction over the nuclear volume. However, the multiple scattering 
of mesons in the nucleus is not accurately accounted for. Our results of the microscopic four-body calculations of the 
r\ c 3 He binding energy are presented in Fig. 2 and compared with the alternative ones in Table II. The calculations 
are performed with the value a = 0.6 GeV for the range parameter of the rj c N potential. 

In summary, we calculated the binding energies of rj c d and rj c 3 He using three- and four-body equations of the AGS 
theory. The results exhibit essential differences to the predictions of the folding model as well as to the variational 
calculations. The characteristic values of the binding energies are about —3 MeV for the rj c d system and amount to 
— 14.5 MeV for a — 0.6 in the 77 c 3 He case. 
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